Weak dynamical localization in periodically kicked cold atomic gases 
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Quantum kicked rotor was recently realized in experiments with cold atomic gases and standing 
optical waves. As predicted, it exhibits dynamical localization in the momentum space. Here we 
consider the weak localization regime concentrating on the Ehrenfest time scale. The later accounts 
for the spread-time of a minimal wavepacket and is proportional to the logarithm of the Planck 
constant. We show that the onset of the dynamical localization is essentially delayed by four 
Ehrenfest times and give quantitative predictions suitable for an experimental verification. 
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Unprecedented degree of control reached in experi- 
ments with ultra-cold atomic gases Q allows to investi- 
gate various fundamental quantum phenomena. A real- 
ization of quantum kicked rotor (QKR) is one such possi- 
bility that recently attracted a lot of attention 
To this end cold atoms are placed in a spatially pe- 
riodic potential Vo cos(2k L x) created by two counter- 
propagated optical beams. The potential is switched on 
periodically for a short time r p -C T, giving a kick to the 
atoms; here T is a period of such kicks. The evolution of 
the atomic momenta distribution may be monitored after 
a certain number of kicks. If the gas is sufficiently dilute 
0, one may model it with the single-particle Hamilto- 
nian, that upon the proper rescaling takes the form |jj of 
the QKR: 
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H = i I 2 + K cos 9 $ ~ n) 
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Here 9 = 2kr,x and time is measured in units of the kick 
period, T . The momentum operator is defined as I = 
ikde, where the dimensionless Planck constant is given 
by k = 8hTk 2 j /(2m). Finally, the classical stochastic 
parameter is K = kVoT p /h . 

The classical kicked rotor is known to have the rich 
and complicated behavior In particular, for suffi- 
ciently large K (^5), it exhibits the chaotic diffusion in 
the space of angular momentum The latter is asso- 
ciated with the diffusive expansion of an initially sharp 
momenta distribution: S(l 2 (t)) = ((l(t) - l(0)f) = 2D cl t 
(dashed line on Fig.QJ. For sufficiently large K, the clas- 
sical diffusion constant may be approximated by K 2 /4 
0. The higher order correction is an oscillatory func- 
tion of the stochastic parameter, i.e., D cl (K) « \K 2 {1 — 
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^h(K) + 2Jl(K)) E3- It was realized a while ago 
|llL UJj that quantum interference destroys the diffu- 
sion in the long time limit and leads to localization: 
S(l 2 (t)) £ 2 at t > t L = D cl /k 2 = £ 2 /D ch where the 
localization length is given by £ = D c \jk. For a large lo- 
calization length £ 3> k, there is a long crossover regime, 
1 < t < it, between the classical diffusion and quantum 




FIG. 1: The momentum dispersion for K — 6.1 and £ = 0.6 
- full line; the classical limit (% — > 0) - dashed line, standard 
weak localization (tE = 0) - dashed-dotted line; the limit 
A2 — * 0, Eq. - long-dashed line. 



localization. 

It was suggested in Ref. ^jj, that the QKR may be 
mapped onto the one-dimensional Anderson localization 
with the long range disorder. The universal long-time be- 
havior of the latter is described by the non-linear sigma- 
model |14| . resulting in the standard weak-localization 
correction [H G3: 5{l 2 {t)) = 2D cl t (1 - 0.75^*7*1) 
(dashed-dotted line on Fig.Q. Notice, that the correc- 
tion is linear in k and non-analytic in time. In an appar- 
ent contradiction with this fact, explicit studies [171 
of the first few kicks show only renormalization of the 
diffusion constant starting from terms ~ k 2 . The aim 
of this paper is to develop a quantitative description of 
the classical to quantum crossover for the QKR that, in 
particular, accounts for these conflicting observations. 

It is known in various contexts that such crossover 
involves an additional time scale, tE, called an Ehrenfest 
(or breaking) time. For a generic quantum mapping, it 
was first shown by Berman and Zaslavski |2Ct \'2\\ , that 
quantum corrections become comparable to the classi- 
cal limit at the time t&. This is the time needed for a 
minimal quantum wavepacket, S9qSIo ~ k, to spread uni- 
formly over the angular direction. Due to the chaotic mo- 
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tion, trajectories diverge as S9(t) — SO^e , where A is the 
classical Lyapunov exponent. For K 3> 1, A = \n(K/2) 
Estimating SIq « KS9 - a typical momenta disper- 
sion after one kick, one finds for the Ehrenfest time: 



(2) 



It is widely believed [HIM mil that this intermedi- 
ate, 1 < < thi time scale is indeed relevant for the 
quantum evolution of classically chaotic systems. This 
observation was put on the quantitative basis in Ref. |23| 
in the context of localization caused by classical scatter- 
ers. In this paper we adopt methods of Ref. [23J to the 
essentially different problem of the QKR. 

In the leading order in k we found for the momentum 
dispersion: 



5(l 2 (t)) =2D cl t- 



3^ 



9{t - U E ) (t - U E ) 



3/2 
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where 9{t) is the step function (long-dashed line on 
Fig. [IJ. At relatively large time, t# <C t < our 
result approaches the standard weak-localization, men- 
tioned above. However, corrections of the order of k 
are absent for t < Me- The delay is caused by the in- 
terference nature of the localization. Indeed, the first 
correction originates from the interference of two closed- 
loop counter-propagating trajectories, see Fig.|2 It takes 
time about t E for classical trajectories passing through 
(almost) the same value of the momentum to diverge and 
take counter-propagating roots. As a result, the inter- 
ference effects are practically absent at smaller times and 
show up only after AtE- 

One may show that the time interval < t < At E is 
protected from higher order weak localization corrections 
as well. For example, in the second order weak localiza- 
tion correction there are two diagrams [24| proportional 
to k 2 9(t — mt E )(t — ratE) 2 with m = 6 and m = 8 corre- 
spondingly. This fact agrees with the perturbative stud- 
ies of the QKR dynamics after a few kicks 0, fill Il2| , 
where no localization effects were seen. (Though the 
classical diffusion coefficient, D c i, is renormalized as an 
analytic function of k 2 .) It is important to mention, 
however, that in reality there is no non-analyticity at 
the point t = 4t E as may seem from Eq. ©. The 
small localization corrections, non-analytic in k, do exist 
for t < 4t E - They are associated with the fluctuations 
of the Ehrenfest time. Since in the quantum mechan- 
ics the minimal separation between the trajectories is 
about 59o ~ y/fc/K, it takes a finite time (~ £e) for 
them to diverge. This time may fluctuate depending 
on initial conditions. The fluctuations are characterized 
by the time scale St E — ^e/^ 2 , where (cf. Eq. (J2J) 
A 2 = (A 2 ) - (A) 2 « 0.82 for sufficiently large K and 
the angular brackets denote averaging over the initial an- 
gle. The interference between rare trajectories, diverging 




FIG. 2: The first quantum correction to the density-density 
correlator: (a) one-loop weak localization diagram; (b) its 
image in the momentum space; (c) semiclassical Hikami box. 



faster than the typical ones, leads to quantum corrections 
at t < 4is of the form: 



5{l 2 {t)) =2D cl t- 



r 



3tt/256 



\6t E ) 3/2 f 



Atp-t 



_ (4) 
where /(0) = 1 and f(x) = 6V2^/T (f ) x -5/2 e -x 2 /ie for 
x 3> 1. Localization correction including effect of the 
Ehrenfest time fluctuations is depicted on Fig. H by a 
full line. It is rather close to the prediction of Eq. © 
(long-dashed line), however the singularity at t — 4Ae is 
rounded. 

Having outlined our main results, we turn now to 
some details of the calculations. One starts from in- 
troducing the exact one period evolution operator as: 

U = exp{i(K/%) cos6>jexp jif 2 /(2£)j. All physical 

quantities may be expressed in terms of the matrix el- 
ements of U n , where n stays for the number of kicks 
(time). We shall be particularly interested in the four- 
point density-density correlator, defined as: 



(5) 



V(l+,l-;l' + ,l'_;uj + ,uj-) 

OO _ 

= Y, (l+\e lu+n e^U n \l' + )(l'_\e^- n 'e J S-U n '\l-) , 

n,n'— 

where \l±) denote momentum eigenstates. First we note 
that averaging over (lo + + u>-) /2 leads to n = n' [lif . 
Then performing the standard Wigner transform, one 
passes to the variables I = (l + + l-)/2 and 9 - the Fourier 
transform of Z+ — Z_ (and similarly for the prime vari- 
ables). Since the most unstable direction of the under- 
lying classical dynamics is along the 9 direction 0, any 
initially smooth distribution quickly relaxes in this direc- 
tion. Thus, averaging over 9 and 9' may be performed. 
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The resulting correlator depends only on the relative mo- 
mentum I — I' and frequency u> = ld + — lu- . Introducing 
finally angle tp as the Fourier image of I — I', one ends 
up with T> — T>(ip;uj). Its classical limit, T> (ip;u)), may 
be found e.g. by using the diagrammatic technique of 
Ref. 0, where in the large K limit it corresponds to the 
family of ladder diagrams. Alternatively, one may show 
from Eq. © that T>o satisfies classical Liouville equation. 
Upon the proper regularization 0, 0, |2j| , that may be 
viewed as a coarse graining in the angular direction, one 
arrives to the classical diffusion propagator: 



V (ip;uj) = (-iu> + D c np 2 



(6) 



This classical limit reflects the diffusion in the momentum 
space: 8(l 2 (t)) = 2D c it, with the diffusion coefficient, 
D c i (K) , studied extensively in the literature [a, 111 LLi| • 

The first quantum correction to Eq. @ is given by 
the one-loop weak-localization diagram, Fig. Et- It de- 
scribes the interference of the two counter-propagating 
trajectories, passing through (almost) the same point in 
the momentum space, Fig. |5Jd. In the Wigner represen- 
tation such correction takes the form: 



dl Q de Q dhdO! r . 
o i "o; t , t/i j 

7T Z7T ^ 



8V{l,6;l,6') -- 

x x(i ,9o- 1 h,9i)[v (i,e-jo,e )Vo(h,9i- 1 i , ,e , )]y (7) 

where u> argument is omitted to shorten notations. The 
operator X stays for the Hikami box [26|. which is given 
by 



X 



exp 



K 2 (69 ) 



■Ai 



Sl o 8d 



IK 



Ak 2 ' " k 

(8) 

where 58q = 9q + 9\ and 81$ = (Iq — li) /2. It is 
clear from this expression that the quantum correction, 
Eq. (7J, is non-zero as long as 5lo58o < k and, there- 
fore, it is proportional to %. The semiclassical Cooperon 
C(Iq, 9o;h,0x) gives the probability of return to (almost) 
the same momentum, l\ w Iq, at (almost) the opposite 
angle, 8\ ss — 9q. If these conditions were strict, such mo- 
tion would be forbidden by the time-reversal symmetry. 
The quantum uncertainty makes it possible. It takes, 
however, a long time to magnify the initially small an- 
gular variation 89q — y/W/K (this estimate as well as 
t5Z — VkK follows directly from Eq. JHJ) up to S9 n w 1, 
when the usual diffusion takes place. 

To take this fact into account we divide the Cooperon 
trajectory onto two parts: the Ehrenfest region, where 
S8 n 1 and the diffusive region with S8 n > 1. We 
denote the corresponding propagators as W and Cq and 
write in the time representation C{t) = j dt' W(t') Co(t — 
2t') (cf. Fig. [2b), where, as we show below, t' ss tE- No- 
tice that the diffusive part of the trajectory is shortened 
by It', leading to C{lo) = W(2uj)C (lu). The diffusive 



Cooperon, Co(lo — h',u), has the same form as Eq. © 
and thus Cq(0;uj) ~ J dip(D c iip 2 — iuj)^ 1 . 

To evaluate propagator W(2ui) in the Ehrenfest 
regime, we define W(z,n) as a probability to reach the 
deviation S9 n = e z during n kicks, starting from an ini- 
tially small variation, S9o ~ y/k/K. According to the 
classical equations (the standard map) 6 n = 9 n -\ + /„ 
and l n — ln-i + K sm8 n -i , the variation evolves as 
S8 n = S9 n - 1 (l + K cos9 n -!) + 2Sl n -!. Since 5l ~ KS9 , 
in the leading order in K 3> 1 the evolution of 89 is given 
by 59 n « 59q n^=o K cos ®i ■ Taking the logarithm of this 
expression, one obtains 




-In \S9 Q 




where the averaging is performed over the initial distri- 
bution of 89q [with the typical scale 89q ~ y/k/K, cf. 
Eq. as well as over dynamics of the fast variable, that 
is 5l n /59 n . For K 3> 1, one may treat cos 6^ after succes- 
sive kicks as independent random variables and employ 
the central limiting theorem to perform the averaging in 
Eq. (JUJ). As a result, 



W(z; n) ~ exp ■ 



(z - ln y/k/K - n\) 2 
2n\ 2 



(10) 



where the Lyapunov exponent A |8|| and its dispersion A2 
are defined as: 

A = (In I if cos 61 1) = \n{K/2) ; (11) 
A 2 = (ln 2 |ifcos6»|)-A 2 = C(3)-ln 2 2«0.82; 

the angular brackets imply integration over 9. The 
Ehrenfest evolution crosses over to the usual diffusion at 
89 n « 1, meaning z w 0. Performing finally the Fourier 
transform as W{lo) = J2 n e lujn W(0, n) and employing the 
definition of the Ehrenfest time, Eq. (J2J), and the fact that 
A2 <C A, one finds: 



W(2w) = exp<^ 2iut E 



2u?\ 2 t E 
X 2 



(12) 



Due to the time-reversal invariance there is an exact 
symmetry between divergence and convergence of the two 
classical trajectories involved in the weak-localization 
correction. This symmetry is illustrated on Fig. [2k . 
Therefore it takes an additional time ~ t E for the two 
distinct semiclassical diffusons to arrive to the point 
la w l\ and #0 « — 6\, bringing, thus, another factor 
VV(2cj). In a slightly different language, one may define 
the Hikami box for a classically chaotic system [23L W\ 
as 2kW 2 (2oj)D c iV 2 , where one factor W(2u>) comes from 
the two legs of the Cooperon, while another originates 
from the two diffusons. Finally, the quantum correc- 
tion, Eq. (7J, reduces to the renormalization of the dif- 
fusion coefficient in the classical propagator, Eq. ©, as 
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D{lo) = D ct + SD(oj) with 

6D(uj)=-^±W 2 (2u) I . \ . . (13) 

Equations (|12fl and l|13|> constitute the main analytical 
results of this work. They describe quantitatively dy- 
namical weak-localization of the QKR with the account 
for the Ehrenfest time phenomena. One may finally ex- 
press the time evolution of the momentum dispersion in 
terms of the frequency-dependent diffusion coefficient. 
The exact relation reads as: 

oo 

— - — % D{u). (14) 

— oo 

Neglecting fluctuations of the Ehrenfest time (A2 — > 
in Eq. \V3 ) and performing frequency and angle inte- 
grations in Eqs. I|14f> . (| 1 31) . one obtains Eq. (|3J) for the 
momenta dispersion. Notice that in this approximation 
the evolution is purely classical at t < AtE- To account 
for the quantum corrections at t < AtE one needs to keep 
the A2 term in Eq. I|12fl . The straightforward integration 
leads to Eq. 10} . 

Our results, Eqs. ®, (@J, are expected to be quantita- 
tively accurate if separations between the relevant time 
scales: 1 < ts < tr, are large enough. This is the case 
when the two dimensionless constants satisfy inequali- 
ties: h < 1 < K . (In the experiments, we are aware 
of [HQ- 7c > 2, and thus ts ^ 1-) Another restriction 
has to do with the dephasing time, j2^. The later 
may ori gin ate from non-perfect periodicity of the kicks 
(noise) |29| |. spontaneous emission as well as from 
many-body collisions between the atoms Whatever 
the nature of the dephasing time, one needs to ensure 
tjg < T0 to observe the weak localization. 
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